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We apply a recently introduced non-perturbative renormalization method to two types of lattice operators: the 
AS = 2 four fermion operator and the heavy-light static axial current, which are relevant for the physics of K 
and B mesons respectively. The results of the non-perturbative calculations of the renormalization constants are 
compared with the corresponding perturbative ones. 



1. Introduction 

We present the results of two applications of a 
general method for non-perturbative (NP) renor- 
malization of generic composite operators [1] (see 
also [2]). The method completely avoids the use 
of lattice perturbation theory (PT), known to 
be a source of uncertainty in the calculation of 
hadronic matrix elements [3]. 

We apply the general method to the AS = 2 
four-fermion operator AS=2 , relevant for Bk, 
and to the heavy-light static axial current A^, 
relevant for fs ■ Both Bk and fs are very impor- 
tant for the study of CP violation. Moreover, the 
AS = 2 operator is important as a preliminary 
study of the operators involved in the AI =1/2 
rule. Indeed, while in the AS = 2 case the chi- 
ral symmetry-breaking Wilson term induces mix- 
ing only with equal dimension operators, com- 
putable in PT, the AI =1/2 case presents mixing 
with lower dimension operators, implying power- 
divergent NP subtractions. 

The NP method introduced in ref. [1] com- 
pletely avoids the use of lattice PT, but not that 
of continuum PT, on which Wilson's OPE is 
based. The NP renormalization condition is ap- 



plied directly to the Green functions of quarks 
and gluons. Given a bare operator 0^ att (a), 
the renormalized operator obtained with the NP 
method, 



latt , 



(1) 



depends on the external states and the gauge, but 
not on method used to regulate the ultra-violet 
divergences. To stress this point, we call the NP 
renormalization scheme Regularization Indepen- 
dent (RI) [4]. The physical operator 



0^{M) = C m {Mln)O m {n) 



(2) 
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is independent of external momenta and gauge 
(up to higher orders in PT and lattice system- 
atic effects) if the Wilson coefficient function 
Cpj(M///) in the RI scheme is calculated with 
the same external momenta and gauge of (1). The 
NP method is valid for any composite operator, 
as long as we can can find a window in the range 
of fj, such that Aqqq <C A* < 0(1/ a) [1]. 

To gain further insight on the size of the higher- 
order corrections, we will compare the NP renor- 
malization constants with one-loop PT, using 
as our expansion parameter both the standard 



bare coupling a 



latt 



g 2 (a)/Air (SPT) and a 



1 



2 2 

H a 


z+ 




Z'2 


Z-i 


0.46 
0.96 
2.47 


0.91(5) 
0.84(3) 
0.85(2) 


0.08(14) 
0.14(7) 
0.22(6) 


0.34(3) 
0.30(2) 
0.33(2) 


0.34(7) 
0.24(4) 
0.23(3) 


SPT 


0.91 


0.12 


0.12 


0.12 


BPT 


0.84 


0.21 


0.21 


0.21 



Table 1 

Values of Z + ,Z\ 2 t 3 for several renormalization 
scales fi 2 a 2 . The PT values are at fi 2 a 2 = 1. 



"boosted" a] 
[3]. 



l.68a\ att , at (] = 6.0 (BPT) 



2. AS = 2 operator 

We consider the renormalization of the four- 
fermion operator 



O 



AS = 2 



75 j 



(3) 



which appears in the weak effective Hamiltonian 
relevant for K°-K° mixing. 

The discretization of the quark action a la 
Wilson, induces a mixing of the operator (3) 
with dimension-six operators of different chirality, 
Osp , Ova, Ospt , determined by CPS symmetry 
(see refs. [5, 6] for their expressions). We define 
the renormalized operator as 



oM 



(aO = 

Z + (O ASz 



(4) 



z x o 



SP 



z 2 o 



VA 



Z a O 



SPT j 



Since there are no AS = 2 operators of dimen- 
sion lower than six, the mixing constants ^1,2,3 
are finite, while the logarithmically divergent Z+ 
renormalizes multiplicatively the subtracted op- 
erator A, ? =2 • The Wilson coefficient function 
sub 

in the RI scheme relating the operator (5) to the 
physical operator, necessary for a NP determina- 
tion of Bk , can be found in ref . [8] . 

We determine the mixing constants Zi with a 
projection method on the four-point amputated 
Green functions [6]. Denoting by Oi, i = 0, . . . , 3 
the operators AS=2 , Osp, Ova, Ospt, we define 
the projectors P 8 - i = 0, . . . , 3 on the amputated 
tree-level Green functions of the operators 



Oi 



Tr P.-A 



(o) 



Sij . The mixing constants 




Figure 1. Chiral behaviour of (O+) as a function 
of X (see text). The NP Z's are at fi 2 a 2 = 0.96. 



Zi are fixed by the condition that the subtracted 
operator A J^ 2 is proportional to the bare free 
operator, i.e. 



TrP fc Af u rV) 



0, k = 1,2,3, 



(5) 



where A A ^~ 2 (pa) is the amputated Green func- 
tion of 0^k =2 ) which is calculated at equal ex- 
ternal momenta p and in the Landau gauge in 
a completely NP fashion from numerical simula- 
tions. Eq. (5) yields a linear non-homogeneous 
system, from which we determine the mixing con- 
stants Zi. The overall renormalization constant 
Z + is then determined by [6] 

Z + ^a)Z- 2 ^a)T^ h = 2 (pa)\ p2= , 



I, 



(6) 



where r^ ) _2 (pa) 



TrP Af u ^(pa), and Z q 



is 

the light-quark renormalization constant, deter- 
mined NP from the conserved vector current [1]. 

We have done a feasibility study of the NP 
method on an ensemble of 36 configurations, on 
a 16 3 x 32 lattice, at (3 = 6.0, and with a single 
value of hopping parameter k = 0.1425 for the 
SW-Clover quark propagator, in the lattice Lan- 
dau gauge. We have performed the calculation at 
different scales fi 2 a 2 , spreading from 0.46 to 2.47. 
In tab. 1, we show the NP values of the Z's for 
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2 2 

H a 


a 




I 


7 


0.46 
0.96 
2.47 


0.022(16) 
0.017(13) 
0.022(14) 


0.23(19) 
0.21(17) 
0.23(17) 


0.78(13) 
0.70(12) 
0.72(12) 


SPT 
BPT 


-0.067(12) 
-0.054(12) 


0.17(15) 
0.17(15) 


0.62(11) 
0.62(11) 



Table 2 



Values of a, [3 and 7 for several renormalization 
scales fi 2 a 2 . The PT values are at fi 2 a 2 = 1. 



several representative values of fi 2 a 2 , and com- 
pare them with the PT result in the same gauge 
and external momenta [6]. Even with low statis- 
tics, for fi 2 a 2 > 1 the Z's have relatively small er- 
rors. As expected the higher-order contributions 
differentiate the mixing constants with respect to 
the PT value, which at one-loop is the same for 
all three [5, 7]. We note that Z2 is less scale- 
dependent than Z\ and Z3, but higher statitics is 
necessary to determine if the variations with fi 2 a 2 
are real or due to statistical fluctuations. 

In order to investigate the effects of the NP 
corrections, we have combined our results with 
the high-statistics computation of the lattice 
matrix elements of the four-fermion operators 
Oi, i = 0,...,3, done in ref. [8]. In fig. 
1, we show the chiral behaviour of (O+) = 
(A°|O AS=2 |A ) latt /(_P 5 ) with the meson at rest 
as a function of X = 8/3/j^m 2 ^ /(-P5) 2 , where 
(P5) = |(0|s75<i|A' }|. Parametrizing the lattice 
matrix element as 

<A'°|0 AS = 2 |A'°> = a + (]m 2 K+ l[p . q) + (?) 

we found that the use of the NP Z's leaves the 
values of [3 and 7 almost unchanged compared to 
those obtained by using the one-loop PT values 
(both the PT and NP values of [3 are compatible 
with zero). In contrast, a changes sign and its ab- 
solute value is reduced by about a factor of three 
in the NP case, becoming compatible with zero. 
This happens for any choice of fi 2 a 2 between 0.46 
and 2.47. The values of a, [3 and 7 for some rep- 
resentative scales of fi 2 a 2 are shown in tab. 2. At 
present, this is the best test of the stability of the 
results which can be done, because the Z's and 
the matrix elements have been computed on dif- 



ferent sets of configurations. Since a should van- 
ish in the continuum limit, we conclude that the 
use of the NP Z's improves the chiral behaviour 
for a large range of values of fi 2 a 2 [6, 8]. 

3. Heavy-light axial current 

Consider now the lattice heavy-light axial cur- 
rent in the static limit 

A v = bj^j 5 q, (8) 

where & is a static heavy quark and q is a prop- 
agating light quark. The NP renormalization of 
(8) is important in improving the accuracy in the 
determination of fs in the static limit, and con- 
sequently also at the physical value of nib- 
We renormalize A v by imposing [9] 

4 att (// a,H^" 1/2 (/"oa)^- 1/2 (H 
^A{pa,Poa)\ po=ll0)P 2 =ll 2 = 1, (9) 

where Zb and Z q are the heavy- and light-quark 
field renormalization constants respectively. Zb 
is defined NP from a conserved scalar current. 
In contrast to the light-light axial current, in the 
static case A v is not finite, so we have no Ward 
Identity to impose [1]. Moreover, we have two 
renormalization scales, /j, and /j,o, instead of one, 
due to the fact that Eichten's static formulation 
of the heavy quark is non-covariant. We will fix 
one of the two scales, e.g. /j,o, and plot Za as a 
function of the other. The physical renormaliza- 
tion constant Za is 

Z A (m b a) = C RI (m6/ / uo,ni6/ / u)^ att ( / u a, / ua),(10) 

where Cpj must be computed, as Z^^ , in the RI 
scheme in continuum PT. The calculation of Crj 
at the next-to-leading order, necessary for a NP 
evaluation of fs , can be found in [9] . 

In fig. 2 we show z\^ obtained with the NP 
method as a function of fi 2 a 2 = p 2 a 2 , with fioa = 
p a = 1,2. They have been calculated on an en- 
semble of 200 configurations, on a 18 3 x 32 lat- 
tice, at [3 = 6.0. The SW-Clover quark prop- 
agator has been computed in the lattice Lan- 
dau gauge, at 3 values of the hopping parame- 
ter, k = 0.1425,0.1432,0.1440, and extrapolated 
to the chiral limit K cr = 0.1455. 
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Figure 2. as a function of [i 1 a' 1 with fi^a = 

1,2. The solid line is from BPT at n^a = 1, while 
the dashed line from BPT at ji^a = 2. 



In order to improve the off-shell matrix el- 
ements, we have also used an improved 0(a) 
heavy-quark propagator, which is simply ob- 
tained by multiplying the unimproved propagator 
in configuration space by [10, 9] 



fit) 



1 



t+i 



(11) 



Moreover, to obtain a better quality of the sig- 
nal for the Green functions, we have made use of 
traslational invariance and inserted the axial cur- 
rent in an arbitrary point, averaging over its posi- 
tions. This method involves the evaluation of the 
heavy-quark propagator between two arbitrary 
times, and the computational effort required is 
reduced by the use of factorization properties of 
the heavy-quark propagator [9]. Although in the 
limit of infinite number of configurations, this 
method is equivalent to the usual method of in- 
serting the current in the origin, we found a much 
better signal. The stability of the signal is shown 
in fig. 2, in which we also report the BPT result, 
calculated in the Landau gauge and with the same 
external momenta of the NP result. We note that 
the NP result lies somewhat below the BPT one, 
over the entire scale range. We would have an 



even greater discrepancy with SPT. This is not 
surprising since the one-loop corrections are quite 
large [11]. 

4. Conclusions 

The results of the NP renormalization applied 
to logarithmically divergent operators, although 
preliminary, are encouraging. In particular, the 
chiral behaviour of the AS = 2 operator seems 
to be significantly improved. This motivate us to 
repeat the calculation, with larger statistics for 
the AS = 2 operator and with unimproved Wil- 
son fermions for the axial current. We plan to 
apply the NP method to power-divergent opera- 
tors, relevant for the Al =1/2 rule. 
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